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Here we consider two phenomena in the vicinity of a black hole deformed by the tidal gravitational 
force of surrounding matter and by a strong magnetic field: equatorial motion of charged particles 
and the decay of a test scalar field. We were able to analyze both phenomena with analytical and 
simple numerical tools, which was unexpected given the low symmetry of the system. We show that 
both the tidal gravitational force and the magnetic field strongly enhance the release of the binding 
energy for the matter spiralling into the black hole. In the presence of the magnetic field, the left 
and right handed rotations of charged particles are not equivalent and for sufficiently large \q\B 
there are stable anti-Larmor orbits very close to the event horizon, although Larmor orbits are only 
stable at some distance from the black hole. The larger the tidal force, the closer the innermost 
stable orbit to the black hole for both types of rotation. It was also shown that the real oscillation 
frequencies of the characteristic quasinormal modes are considerably suppressed by the tidal force. 
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I. INTRODUCTION 

Although the interior of black holes is invisible, there 
are reports on observations of astrophysical black holes 
owing to their intense interactions with the surrounding 
matter, which are manifest around galactic black holes 
interacting with the Active Galactic Nucleus. Accretion 
disks create gravitational tidal forces as well as electric 
and magnetic fields around black holes. Accurate physi- 
cal descriptions of such processes usually involve complex 
numerical modeling using supercomputers and present a 
number of challenges connected to explanations of vari- 
ous observed processes, including for example, the origin 
of jets around black holes. Here we adopt a more the- 
oretical framework, which, although far from being an 
accurate astrophysical description, could be a bridge be- 
tween the mathematical theory of black holes and the 
numerical modeling in astrophysics. 

An important constituent of a black hole's environment 
is a magnetic field, which can be as strong as 10 4 — 10 8 G 
near the black hole [l| . The magnetic fields near a black 
hole can extract its rotational energy (Blandford-Znajek 
effect 0), affect the quasinormal spectrum Q and Hawk- 
ing radiation Q, jik and induce an external electric 
charge on its surface 3|. At the same time, there is no ev- 
idence that the magnetic field could be strong enough to 
deform the geometry of a black hole significantly. Thus, 
the most interesting effects associated with the presence 
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of the magnetic field occur on charged particles and fields. 
Indeed, the factor which stipulates the Lorentz force dur- 
ing the particle's motion is qB//i [3], where fi is the mass 
of the particle and B is the asymptotic value of the mag- 
netic field strength. 

Another essential factor representing the black hole en- 
vironment is a tidal gravitational force produced by all 
the surrounding matter. Unlike magnetic fields, tidal 
forces can be strong enough to deform (albeit mini- 
mumly) the geometry of a galactic black hole as well as 
black holes of stellar masses having a star companion. 

The perturbative solution to the Einstein-Maxwell 
equations which allows us to study, the magnetic field 
and the tidal gravitational force, was found in [5[ by Pre- 
ston and Poisson with the help of the light-cone gauge 
formalism |6]. They considered a system consisting of a 
black hole and a mechanical structure (a giant torus or 
long solenoid) at some distance around it, which was the 
source of the gravitational tidal force and the strong mag- 
netic field. Thus, the background which we shall study 
here, the Preston-Poisson space-time, has three parame- 
ters: mass M, magnetic field B and tidal force £. 

We investigated two basic properties of particles and 
fields in the vicinity of such a magnetically and tidally 
deformed black hole: the proper oscillation frequencies 
in the black hole's response to the perturbation, termed 
quasinormal modes, and the motion of the charged parti- 
cles. The quasinormal modes do not depend on the kind 
of perturbation but only on the parameters of space-time 
and are therefore termed the fingerprints of the space- 
time. The quasinormal modes are studied in various dis- 
ciplines such as black hole physics, gauge/gravity corre- 
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spondence, gravitational wave astronomy Q, cosmology 
Q and quantum gravity @- 

Particle motion and the quasinormal modes have been 
studied for black holes immersed in a strong magnetic 
field i, 1, 1, 0, [Ht [3, 0, 0, pl.We 
previoulsy investigated motion of neutral particles in the 
vicinity of the Preston-Poisson space-time [I?} • Further 
studies on particle motion, accretion and thermodynamic 
properties around black holes in a magnetic field was per- 
formed in pH, [H, [H, HI- In the presented work 
we estimate quaisnormal modes of a test scalar field prop- 
agating in the background of the Preston-Poisson black 
hole [5j, when the magnetic component of the deforma- 
tion vanishes. This is a good approximation, since the 
gravitational tidal force is usually much stronger than 
the deformation caused by the magnetic field. In the 
ekional regime, quasinormal modes usually have a kind 
of universal behavior for different boson fields, such that 
the spectrum of the neutral scalar field considered here 
could reveal some features of the gravitational spectrum. 
The latter is much more difficult for analysis due to the 
inseparability of variables in the perturbation equations. 

The motion of charged particles is essentially influ- 
enced by the magnetic component. Therefore, in the first 
part of our work, when considering equatorial motions of 
charged particles, we have taken both the gravitational 
tidal force and the space-time deformation due to the 
magnetic field of the torus into consideration. The an- 
gular variables in the Hamilton- Jacobi equation in the 
Preston-Poisson space-time can be decoupled only in the 
equatorial plane, by which our analysis of particle motion 
is limited. In particular, we show that the tidal gravita- 
tional force as well as the magnetic field strongly enhance 
the release of the binding energy for a particle spiralling 
into the black hole. The region of stability is also sig- 
nificantly affected by the tidal force and magnetic field 
and is qualitatively different for left and right handed ro- 
tations. We shall discuss these issues in detail in Sec. 
IV. 

The scalar field equation also does not allow for a com- 
plete separation of variables in the general case. Never- 
theless, we used the fact that the astrophysically moti- 
vated low-laying quasinormal modes are " localized" near 
the peak of the effective potential and therefore "aver- 
aged" the tidal force near the peak by its value in the 
maximum. This allowed us to decouple variables and es- 
timate the dominant quasinormal modes for the first and 
higher multipoles £ > 1 . We have shown that the real os- 
cillation frequencies are considerably suppressed by the 
tidal force in the vicinity of the magnetic field. 

The paper is organized as follows. In Sec. II we 
summarize some basic relations for the Preston-Poisson 
space-time. In Sec. Ill the effective potential for 
the equatorial motion is deduced on the basis of the 
Hamilton- Jacobi formalism. Properties of particle's mo- 
tion are studied in Sec. IV with consideration of various 
particular cases, such as zero tidal force, zero magnetic 
field, and the approximation of the neglected "geometric" 



influence of the magnetic field. Sec. V is devoted to the 
decoupling of variables in the test scalar field equation 
for the Preston-Poisson metric with vanishing magnetic 
field. In Sec. VI, the angular part of the scalar field 
equation is studied. Sec. VII presents calculations of 
quasinormal modes for the above case. In Sec. VIII we 
summarize the obtained results and discuss the questions 
which remain for further study. 



II. PRESTON-POISSON METRIC 

The Preston-Poisson space-time describes the system 
consisting of a large mechanical structure, such as a gi- 
ant solenoid or a torus which surrounds a black hole 
and produces an asymptotically uniform magnetic field 
of strength B. The mass of the structure is M' and the 
radius is of order ~ a, while the mass of the black hole 
is M . The electromagnetic four- vector has the following 
form, 
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(1) 



where <f>^ — (0,0,0,1). Further, it is implied that the 
perturbation created by the magnetic field is small, i.e. 



r 2 B 2 < 1, 



(2) 



where r is the distance from the black hole, and only the 
interior of the mechanical structure is under considera- 
tion r < a. For our purposes it is sufficient to study the 
inner region of the system which starts at the black hole 
horizon (r^ = 2M) and ends far from the black hole, still 
being far from the edges of the torus, 



ri<r«a. 



(3) 



The latter condition can always be fulfilled, because the 
torus is supposed to be situated in the region of the weak 
gravitational field of the black hole, 



M 



< 1. 



(4) 



As r 2 B 2 <C 1 and r < a, it is implied that a 2 B 2 <C 
1, though the relative scales of M/a and a 2 B 2 can be 
arbitrary. First of all, we are interested in the case 



M/a < a 2 B 2 , 



(5) 



i.e. in the situation when there is an asymptotic region 
M <C r < a, where the influence of the magnetic field 
on the space-time geometry is negligible. The mechani- 
cal structure of mass M' produces the gravitational tidal 
force, parameterized by £ , near the black hole, 



£ 



M' 



(6) 



In the above approach the tidal force can be much larger, 
of the same order or much smaller than B 2 . 
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FIG. 1: a) Left panel: Effective potentials U^ ff for L = 3, fx = 0.1, £ = 0.001, B = 0.01, M = 1, q = -10, 0, +10 (from bottom 
to top), b) middle panel: Effective potentials for L = 3, fx = 0.1, £ = 0.001, 0.005, 0.01 (from bottom to top), B = 0.01, M = 1, 
q = +10, c) right panel: Effective potentials for L = 3, fx = 0.1, £ = 0.001, 0.005, 0.01 (from bottom to top), B = 0.01, M = 1, 
(7 = -10 
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Preston and Poisson used the light-cone gauge for con- 
structing the perturbed metric, which is adapted to in- 
coming light cones v = constant that converge toward 
r = 0. For zero tidal force and magnetic field, v takes its 
Schwarzschild value v = t + r + 2Mln(r/2M - 1). In the 
(v, r, 9, coordinates the Preston-Poisson metric has 
the following form 

ds 2 = g vv (r, 6)dv 2 + 2dvdr + g v $(r, 9)dvd9+ 



FIG. 2: Effective potential Ul ss for L = 0,0.1,0.4,0.8 (from 
top to bottom on the right side), M = 1, fx = 0.01, B = 0.01, 
£ = 0.001, q = -5. 



gee(r,9)d9 2 + 5<w ,(r, 6>)d0 2 



(7) 



where 



J 



5ot = -/ - \B 2 r{3r- 8M) - i.B 2 (3r 2 - UMr + 18Af 2 )(3 cos 2 6> - 1) + £ (r ~ 2M) 2 (3 cos 2 9 - 1) + 0[S 4 , £ 2 ], (8) 
9 9 



7r« 



= l + 0[B\£\ 



= -B 2 r 2 (r - 3M) sin (9 cos (9 - 2£r 2 (r - 2M) sin 6» cos 6» + (9[S 4 , £ 2 ] , 



(9) 
(10) 



r 2 - -B 2 r 4 + is 2 r 4 (3cos 2 - 1) + B 2 M 2 r 2 sin 2 9 + £r 2 (r 2 - 2Af 2 ) sin 2 9 + 0[S 4 ,£ 2 ], (11) 
9 9 



g vip = r 2 sin 2 9 - ^B 2 r 4 sin 2 9 + ^B 2 r 4 sin 2 0(3 cos 2 - 1) - B 2 M 2 r 2 sin 4 - £r 2 (r 2 - 2M 2 ) sin 4 + 0[B 4 , £ 2 ], (12) 
9 9 



and/(r) = 1— (2M/r). The above metric was obtained in 
5] by the perturbation of the Enstein-Maxwell equation 
in orders of (£,B 2 ). The parameter £ is the Weyl curva- 
ture, that is, the tidal gravitational field, of the asymp- 
totic space-time measured by an observer co-moving with 
the black hole in the region M < r < 1/5. The per- 



turbed event horizon is given by 

r h = 2M(1 + -M 2 B 2 sin 2 i 
3 



(13) 



It is essential that the event horizon is affected by B and 
not by £. 
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Although in the above space-time only the dominant 
order of the magnetic field is considered, the same rela- 
tion as for the Ernst black holes takes place 



B 



91 M 

10 ^ 



(14) 



where M and M© are the mass of a black hole and of the 
sun respectively and Bq is the external magnetic field 
in units of gauss. From the above relation one can see 
that the magnetic field deforming the space-time geom- 
etry significantly would be as strong as B — 10 12 G for 
galactic black holes with mass M ~ 10 9 M Q . 



and E and L are constants of motion, the particle's en- 
ergy and angular momentum respectively (po = —E, and 
P3 = L) while pi = S r and pi — Se are functions of r and 
9 respectively. In the general case the angular variables 
in the Hamilton- Jacobi equations cannot be decoupled. 
Nevertheless, one can decouple variables for motion in 
the purely equatorial plane 9 = ir/2. 

In the equatorial plane the Hamilton- Jacobi equation 
for the metric (fl6l) takes the form 



9vv 



'ds\ 2 
K dt) 



'dS\ 2 
y dr ) 



III. EQUATIONS OF MOTION AND THE 
EFFECTIVE POTENTIAL 



The transformations 



v = t + F{r), F'(r)g vv (r, 6) + 1 = 
reduce the metric ([7]) to the following form 



(15) 



ds 2 = -g vv (r, 9)- 1 dr 2 + g vv (r, 6)dt 2 + g u {r, 9)d<p 2 + C, 

(16) 

where C is the part of the metric which vanishes on the 
equatorial plane 9 = 7r/2 as it contains d0. The metric 
coefficients g vv and g^ depend on the radial and angular 
coordinates r and 9, the black hole mass M, magnetic 
field B and the tidal force parameter £ . The above metric 
is diagonal on the equatorial plane and convenient for 
further analysis of geodesic motions with the help of the 
Hamilton- Jacobi formalism. 

The Hamilton- Jacobi equations for the time-like and 
null geodesies in curved space-time are 



1 



'""I—- qA, 



as 

dx v 



dS_ 
~ds' 



(17) 



Here s is an invariant affine parameter. The action is 
S = -^ 2 s -Et + Lcb + S r (r) + S e (9), 



OS 

- rg4 A l ~ 2( i A ml gi = ■ (is) 

Implying the normalization r = (is, where r is the 
proper time, the first integrals of motion are 



dr 
ds 



L + ^Bqg^ 



9<P<P 



(19) 



The qualitative analysis of the motion can be made by 
considering the effective potential 



jjZ 2 



(20) 



The effective potential for various parameters is shown 
in fig. [TJ Let us take M = 1 and consider the situation 
when, L 3> Bq and at the same time B <C 1, £ <C 1, 
fi <C 1. Expansion in terms of small fi, B and \q\B / L 
shows that for sufficiently large values of L, the effective 
potential has a local minimum at 



2L 



1/2 



r 



WqB 



3)£- 



(21) 



For small values of the angular momentum L the effec- 
tive potential is a monotonically growing function of r. 
For more general cases, which are not limited by small 
values of /i, it may be useful to introduce the angular 
momentum, charge and energy "normalized" by mass: 

l^l/h, q—tq/fJi, E—^-E/jj,. 



Then, the effective potential can be written as 



v eff 



U 2 

A* 2 



-9v 



1 



(22) 



Further, we shall use this form of the effective potential 
and the " normalized" values of the energy, angular mo- 
mentum and charge. 
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The effective potential, and consequently the motion of 
particles, is not the same for opposite charges of particles 
<7, that is, the left-hand and right-hand rotations are not 
equivalent due to the presence of the magnetic field: the 
rotation which corresponds to the Lorentz force directed 
from the black hole is called the anti-Larmor rotation, 
while the rotation with a Lorentz force directed towards 
the black hole is the Larmor one. In the next section, we 
shall show that negative (positive) charge and the pos- 
itive (negative) angular momentum corresponds to the 
Lorentz force directed from the black hole (anti-Larmor 
motion), and, when the signs of the charge and of the 
momentum coincide, the motion is Larmor. 

From fig. [1] one can see that the effective potential 
of positively charged particles is higher than that of the 
negatively charged ones at all r. Thus, the magnetic 
field allows negative particles to penetrate the barrier at 
energies smaller than those for the Schwarzschild black 
hole. 

Now, we are in position to study properties of the equa- 
torial motion of charged particles. 



IV. PROPERTIES OF PARTICLE MOTION 



B = in formulas ([7]) for g vv and and keep B ^ 
in the generalized derivatives of the Hamilton- Jacobi 
equations (fl~7|) . This approximation is quite good, taking 
into account that the magnetic field around the black hole 
decays quite quickly with r, so that its overall influence 
on the geometry can indeed can be neglected for small 
qB. 

The parameters of the circular orbits can be deter- 
mined from the requirements: 



dr 
ds 



d 2 r 

^ = ' 5? 



o. 



(23) 



These equations mean that on the circular orbits one has 



dVeff 

V eff = E, -£1 = 0, 



(24) 



Here, we shall concentrate on circular orbits of charged, 
massive particles. The case of massless and neutral par- 
ticles was considered, though in a different coordinate 
system, in (T^. 



that is, circular orbits take place only on the turning 
points of the effective potential. Next, let us consider 
the two particular cases: 1) vanishing tidal force and 2) 
vanishing magnetic field. 

1. The case £ = 



A. The case where the influence of the magnetic 
field on the black hole's geometry vanishes 

We shall first neglect the deformation of the black hole 
geometry due to the magnetic field, that is, we shall take 



In the case of vanishing tidal force, £ = 0. Solutions to 
the system of equations (|24|) allow us to find the values 
of energy E and momentum L as a function of the rest 
of the parameters B, M, q and coordinate r, 



J 



E z 



(r - 2M) 2 [B 2 q 2 {r - 2M)r 2 + 2r - 6M T Bq^Jr 2 {B 2 q 2 r A + AM 2 {B 2 q 2 r 2 - 3) + M (4r - 4B 2 q 2 r 3 )) 

2r(r - 3M) 2 



(25) 



L 



± v / r 2 (B 2 q 2 r i + 4M 2 (B 2 q 2 r 2 - 3) + M (4r - AB 2 q 2 r 3 )) - BMqr 2 

6M - 2r 

I 



(26) 



When B = 0, the above expressions ([25] 126 ]) for the en- 
ergy and momentum reduce to their Schwarzschild values 



E 



(r - 2M) 2 



Mr 



r(r-3M)' 1 ^r - 3Af ' (2?) 

Note that L, E and q are "normalized" by mass /it. 
Under positive g, the upper sign corresponds to the anti- 
Larmor motion (A), while the lower sign, to the Lar- 
mor one (L). Negative q gives opposite correspondence 



for Larmor and anti-Larmor rotations. Indeed, asymp- 
totic behavior at large (far from the black hole) r, gives 
the following asymptotic values for the energy and mo- 
mentum: 



E\ -> 1, L A -> -\Bqr 2 , 



El B 2 (3M 2 + r 2 ) q 2 + 1, 



r»M, 



r > M, 



(28) 



(29) 
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L L ^ -Bq (6M 2 + 2rM 



r > M. (30) 



The Larmor motion is a kind of cyclotron rotation in 
a uniform magnetic field, when the magnetic field is di- 
rected to the center of the orbit. The Larmor motion 
does not require a black hole, so that a black hole just 
perturbs the existing Larmor orbits. On the contrary 
to the Larmor rotation, anti-Larmor motion takes place 
only in the presence of the black hole, when the Lorentz 
force is directed outwards from the center. 



2. The case B = 

The case of vanishing magnetic field is more trivial as 
the motion of charged particles is qualitatively the same 
as those of the neutral ones. The energy and momentum 
expanded until the first order in £ have the following 
form: 

e2 £ (8M 2 - 7rM + 2r 2 ) (r - 2M ) 2 



(r - 3M) 2 



(r - 2M) 2 
r(r - 3M) 



0(£ 2 



(31) 



L = ±- 



£ (-6M 4 - 2rM 3 + 9r 2 M 2 - 5r 3 M + r 4 1 r 



2\[M(r - 3M) 3 / 2 



+ 0{£ 2 ). 



Mr 



Vr - 3M 



-± 



(32) 



It is evident that the above expressions reduce to their 
Schwarzschild values as £ = 0. 



3. The case B ^ and £ / 

For a more general case B^O and £ ^ (but still in 
the approximation of the neglected geometrical influence 
of the magnetic field onto the metric), the explicit ex- 
pressions for the energy and momentum are rather cum- 
bersome, therefore figs. (JSJIH]) were used to depict the 
dependence of energy and momentum on the radius of 
circular orbits. 

In figs. <j3j> and ([5]) one can see that the energy E and 
absolute values of the momentum L increase when the 
tidal force £ increases. Usually, the effective potential 
has two extrema: one at the maximum of the potential 
barrier (corresponding to an unstable orbit and not far 
from the black hole) and one at the minimum which is 
further from the black hole (corresponding to a stable 
orbit). At some minimal values of the angular momen- 
tum, these two extrema coincide, which corresponds to 
the innermost stable circular orbit. Thus, the minimum 



of the absolute value of angular momentum L in figs. 

is located at the innermost stable circular orbit. 
As it is shown in figs. (j3]) and ([5]), the larger the tidal 
force £ , the closer the innermost stable orbit to the black 
hole. This could be explained by the fact that the tidal 
force £ of the surrounding structure is opposite to the 
gravitational attraction of the black hole, such that the 
"effective" gravitational attraction becomes weaker and 
particles can come closer to the black hole yet still remain 
in a stable orbit. 

The dependence on charge q is more complicated as the 
situation depends on the direction of the Lorentz force. 
Thus, the anti-Larmor motion (when the Lorentz force 
is directed outward from the black hole) corresponds to 
positive q and negative L or to negative q and positive 
L. The energy and momentum for various q are shown in 
figs. (|4]) and (j6|). The positive and negative angular mo- 
menta show opposite dependence of the energy on charge 
q: for negative L, lines of energy for negative q always lay 
above those for positive q, while for positive L the situ- 
ation is reversed, that is, negative charges correspond to 
smaller E than positive ones. The same dependence takes 
place for the absolute value of the angular momenta. In 
summary, the energy E and the absolute value of the an- 
gular momentum \ L \ grow with the increasing absolute 
value of charge \ q \ for Larmor orbits and decreases for 
anti-Larmor ones. As the magnetic field B is coupled to 
the charge q, the dependence on B is qualitatively sim- 
ilar to the dependence on q. As the dependence on the 
magnetic field for neutral particles was considered earlier 
in , we focused on the dependence of various quanti- 
ties on the charge and the tidal force. For the relatively 
small q (when qB is still small) and £ , considered in figs. 
([3]) and © and in Tables I and II, one cannot tell the 
difference in the location of the innermost stable Larmor 
and anti-Larmor orbits. That is why we have only one 
set of values of r in Tables I and II: the corrections to 
those values are of the order 0(B 4 , £ 2 , B 2 £). This is 
certainly not the case for large q for which the location 
of Larmor and anti-Larmor orbits are quite different. We 
shall consider the regime of large q in the next subsection. 

In order to learn how much energy could be extracted 
from a particle slowly spiralling toward the center, it is 
necessary to know what the energy of the particle at the 
innermost stable circular orbit is. For this purpose it is 
convenient to use a function called the binding energy, 
which is the amount of energy released by the particle 
going from the stable circular orbit located at r to the 
innermost stable circular orbit r^. Thus, the binding 
energy (in percent) is 



Binding Energy = 100 



E(r)~E(r ls ) 
E(r) ■ 



(33) 



The binding energy can tell us how much energy the 
matter (for instance an accretion disk) will release before 
plunging into the black hole. For the Schwarzschild black 
hole (taking M = 1) the energy on the innermost stable 
orbit is E(ri S ) = \/8/9, so that for a particle coming 
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q 


r 


L+ 


L- 


E- 


E+ 


-10 


5.2745 


3.64979 


-3.78766 


0.99508 


0.95882 


-3 


5.2841 


3.69291 


-3.73427 


0.98163 


0.97076 


-1 


5.2849 


3.70624 


-3.72003 


0.97792 


0.97432 





5.2850 


3.71308 


-3.71308 


0.97612 


0.97612 


1 


5.2849 


3.72003 


-3.70624 


0.97432 


0.97794 


3 


5.2841 


3.73427 


-3.69291 


0.97432 


0.97794 


10 


5.2753 


3.78766 


-3.64979 


0.95882 


0.99508 



TABLE I: Energies E + and _E_ and the angular momenta 
L+ and L- at the innermost stable circular orbit for various 
values of charge q; B = £ — 10~ 3 , M = 1 



q 


r+ 


r_ 


L+ 


L- 


E- 


E+ 


-10 


5.27452 


5.27527 


3.64981 


-3.78768 


0.995088 


0.95882 


-50 


5.04704 


5.11032 


3.49905 


-4.19065 


1.08802 


0.90389 


-100 


4.86491 


4.61837 


3.46307 


-4.87639 


1.23937 


0.85679 


-1000 


2.74943 


4.32074 


5.20535 


-24.8736 


5.90001 


0.59092 



TABLE III: Energies E + and £_ and the angular momenta 
Z/+ and L- at the innermost stable circular orbit for various 
large values of charge q; B — £ — 10 -3 , M = 1 



The general case with all non-vanishing 
parameters 



£/Hn 3 


r 


L+ 


L- 


E- 


E+ 





5.9997 


3.47018 


-3.45818 


0.94119 


0.94445 


0.5 


5.5277 


3.60380 


-3.59086 


0.95811 


0.96156 


1 


5.2849 


3.72003 


-3.70624 


0.97432 


0.97794 


1.5 


5.1213 


3.82919 


-3.81458 


0.99046 


0.99425 


2 


4.9982 


3.93531 


-3.91988 


1.00684 


1.01080 


2.5 


4.8994 


4.04064 


-4.02438 


1.02366 


1.02779 


3 


4.8168 


4.14672 


-4.1296 


1.04105 


1.04537 



TABLE II: Energies E+ and E- and the angular momenta 
L+ and L_ at the innermost stable circular orbit for various 
values of charge £; B = 10~ 3 , q = 1, M = 1 



from the asymptotic region (E = 1), the binding energy 
is about 5.7 percent. For the Preston-Poisson space-time 
there is no asymptotically flat region, so that, the binding 
energy released at the transition from a given stable orbit 
r (instead of infinity for the Schwarzschild case) to the 
innermost stable orbit r^ s is meaningful. Figs. ([7]) and 
(|8]) depict how the binding energy depends on the charge 
q, tidal force £ and magnetic field B as a function of r. 

In fig. (J7J) the line for negative L lay above the line 
for positive L, if q is negative, and below, if q is positive. 
Thus, the Larmor orbits have greater binding energy than 
the corresponding anti-Larmor ones. Comparison with 
the blue curve for the Schwarzschild case (B = £ = 0) 
in fig. (O demonstrates that the magnetic field and its 
coupling with the charge strongly enhance the release of 
the binding energy. This property was also observed in 
13] for the black hole in the absence of the tidal force. 
The tidal force also strongly increases the binding energy, 
as can be seen from comparison with the Schwarzschild 
case in fig. ((SJ). There, the negative momenta correspond 
to a slightly smaller binding energy than the positive one 
in fig. ©, where the small difference is due to the small 
Bq. Thus, the Larmor and anti-Larmor orbits are almost 
indistinguishable in the regime of tiny Bq. In the next 
subsection we discuss the regime of large charge, when 
Bq is not necessarily small. 



When the charge q is large enough, so that the term 
Bq is no longer small, one cannot safely neglect the geo- 
metric influence of the magnetic field, because even sub- 
dominant terms containing Bq can be quite large in com- 
parison, for instance, with the tidal force corrections ~ £. 
Thus, we considered the full metric ([7]) with a nonzero 
magnetic field both in the metric coefficients and in the 
Hamilton- Jacobi equation for large q. As shown in Table 
III, the locations of Larmor and anti-Larmor orbits al- 
ready differ a lot for Bq ~ 1. At the same time, the data 
for Bq ~ 10~ 2 (the first line in the Table III) still almost 
coincide with the approximation of the neglected "ge- 
ometric contribution" of B considered above. Another 
distinction from the above small charge approximation 
is that, while for small Bq the energy at the innermost 
stable orbit is only slightly smaller for the anti-Larmor 
orbits than for the Larmor ones (Tables I and II), for 
moderate Bq, this difference is very large (see last line 
in Table III) and the Larmor orbits may have one order 
larger energy and momentum than the anti-Larmor ones. 

In the case of large q, the binding energy of Larmor or- 
bits is not always greater than the anti-Larmor ones as it 
depends on a stable orbit from which the particle plunges 
into the black hole. For more distant stable orbits (see 
fig. ©) the binding energy of the anti-Larmor orbits can 
be larger than the corresponding Larmor ones. If one re- 
members that the Larmor motion is a cyclotron rotation 
in a uniform magnetic field perturbed by a black hole, 
then, it becomes clear that in the case of large Bq the 
"perturbation" by the black hole is relatively small, so 
that the rotation is almost purely cyclotronic, especially 
at a relatively large distance from the black hole. 

Finally, let us discuss the radial stability of the consid- 
ered circular orbits. The condition of the radial stability 
can be written as: 

d 2 U rf 

- , > - stability condition. (34) 
or 

From numerical study of the above inequality (when U e f / 
is given by its most general expression (|22p). we conclude 
that at sufficiently large \q\B, the region of stability of 
anti-Larmor orbits can approach the event horizon, 

r > rh - anti-Larmor, high \q\B, (35) 
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FIG. 3: Energy E- (left) and momentum L_ (right) as a function of circular orbit radius for various values of £ — 
0, 1/300, 1/200, 1/100 (from bottom to top for energy E- and from top to bottom for the negative angular momenta L_), 
B = 1/100, M = 1, q = 1. 




FIG. 4: Energy E- (left) and momentum L_ (right) as a function of circular orbit radius for various values of the particle 
charge q = -5(green), 0(red), 5 (blue) £ = 1/100, B = 1/100, M = 1. 





FIG. 6: Energy _B+ (left) and momentum L+ (right) as a function of circular orbit radius for various values of the particle 
charge q = -5(green), 0(red), 5 (blue), £ = 1/100, B = 1/100, M = 1. 
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FIG. 7: The binding energy as a function of r for q = -10 (left), q = (middle) and q = 10 (right), 8 = B = 10" 3 , M = 1; red 
is for negative L and green is for positive L, blue line is for the Schwarzschild orbit. 



Binding Energy (%) 



nding Energy (%) 




Binding Energy (%) 



FIG. 8: The binding energy as a function of r for q = 1, £ = 0.5 x 1CT 3 (left), £ = 1.5 X 1CT 3 (middle) and £ = 3 x 1CT 3 
(right), B = 10~ 3 , M = 1; red is for negative L and green is for positive L, blue line is for the Schwarzschild orbit. 



Binding Energy (%) 




FIG. 9: The binding energy as a function of r for q = —1000, 
£ = B — 10~ 3 , M = 1, red line corresponds to the orbit with 
a negative momentum (Larmor for q < 0), green corresponds 
to a positive (anti-Larmor) 



while Larmor orbits are only stable for r 2; 4.2M for zero 
tidal force £ = 0, approaching r gl 3.6M at £ = 1/50M, 

r I 4.2M(£ 0) \ Larmor, high Igl-B. (36) 
~ ^ 3.6M(£ = 1/50M) y 

This is in good agreement with the r ^ 4.3M limit 
observed in [3| for the Ernst solution in the absence of 
the tidal force. The small difference between ((551) and Q 
is apparently due to the neglected contribution of orders 
B 3 and higher in the Preston-Poisson metric. The new 
feature here is that the tidal force can considerably ex- 
pand the region of stability of the Larmor orbits making 
them closer to the black hole horizon. This can also be 
seen from the Table II. 



V. SCALAR FIELD PERTURBATIONS IN THE 
PRESTON-POISSON SPACE-TIME 



The aim of our study of a test scalar field in the 
Preston-Poisson space-time is to find the proper oscilla- 
tion frequencies (quasinormal modes) which dominate in 
the response to the perturbation at late time. As shown 
in [23j . the quasinormal modes of black holes, when con- 
sidered in the astrophysical context, are essentially inde- 
pendent of the behavior of the master wave equation far 
from the black hole. The explanation is straightforward: 
quasinormal modes are poles of the reflection coefficient 
of the scattering process which occurs near the peak of 
the effective potential, so that the low-laying, dominating 
modes are " localized" near the maximum of the effective 
potential. For the Schwarzschild solution this happens 
at r = 3M, while for the Preston-Poisson space-time this 
value is only slightly corrected. The wave equation does 
not have physical meaning at r = oo, because for as- 
trophysical processes of this kind, "infinity" effectively is 
situated at the distance which is much larger than the 
radius of the black hole. In our case, "infinity" is at least 
quite a few times larger than the black hole radius and 
still far from the edge of the torus, <C r -€1 a. We 
have considered the same spatial region when analyzing 
particle's motion. 

For a test neutral scalar field, influence of the magnetic 
field B can usually be neglected, because the energy den- 
sity of the magnetic field is much smaller than that of the 
gravitational one. When B = 0, the coefficients of the 
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Preston-Poissin metric ([7]) are simplified as follows: 

9w = -f + S(r - 2M) 2 (3cos 2 -1), 
Qvr — 1? 

g v g = -2£r 2 (r - 2M) sin 9 cos 9, 
gee = r 2 + £r 2 (r 2 - 2M 2 ) sin 2 9, 



; sin 2 9 - £r 2 {r 2 - 2M 2 ) sin 4 



(37) 
(38) 
(39) 
(40) 
(41) 



The Klein-Gordon equation in a curved space-time has 
the following form 



-9 



.9* 



= 0. 



In the general case, the variables in the Klein- Gordon 
equation cannot be decoupled for the above metric. How- 
ever, a further assumption can remedy the situation. We 
shall assume, that as the torus is situated far from the 
black hole M <C a, the tidal force £ , being relatively 
small, acts almost homogeneously in a small region where 
the low laying quasinormal modes are "localized", that 
is near the peak of the effective potential 



: PM, p 



3' 



0. 



(42) 



P = 



y/9 + 1(1 + 1)(14 + 91(1 + lj) - 3 
21(1 + 1) 



In other words, the tidal force at the potential peak does 
not change much when one slightly moves away from the 
peak r s» r max . Thus, we assume that M -C a holds 
"with a margin", being also r max < a. In this case, be- 
ing interested only in the region near the peak, we can 
expand the perturbed part of the Preston-Poisson met- 
ric (I37H41P (i.e. the terms containing £) in powers of 
(£, (r — r max )). Further we neglect the corrections to the 
Schwarzschild metric of order £ (r — r max ) and higher, 
which means that, as we assumed, the tidal force barely 
changes near the peak. Although such an approximation 
does not allow us to make accurate calculations of quasi- 
normal spectrum, we can still estimate the dominating 
modes numerically. In reality the tidal force is smaller 
on the left (at the black hole side) of the peak and larger 
on the right of it, so that the used expansion in powers of 
(£ , (r — r max )) is a kind of "averaging" of the tidal force 
in the small region of the localization of the dominant 
modes. 

In the above approximation the metric can be reduced 
to the diagonal form 

ds 2 = g vv dt 2 —dr 2 + geedO 2 + g w d(p 2 , (43) 

Qvv 

through the following transformations, 

v = t + F(r), F'(r)g vv (r,9) + 1 = 0, (44) 



where 

9w = ~f, 9ve = 0. (45) 

The determinant of the diagonal metric (|43|) is 

V-9 = V9ee9<ptp ~ r 2 sin6» + 0(£ 2 ). 

Due to the Killing vectors in t and cj>, the stationary 
ansatz for the perturbation is implied: 

9(t,r,0,(p) =e vmv - lut ${r,9). 

The wave equation can be written as follows 



or \ \ r J or 



S$ = 0, 



where 



5$ = 



Expansion of ./^- and * &££. i n powers of (£Ar 
V Sff y gee * \ i \ 

)) yields 



^ = sin 9-{p 2 - 2)£M 2 sin 3 9 + 0(£(r~r ' c ' 2 



^max ) ; & ) i 

gee 

— = -^ + (p 2 -2)M 2 £sm9 + 0(£(r-r max ),£ 2 ). 
g vv sin 9 

Then, the angular part can be rewritten as 
1 d ( . M\ (p 2 -2)£M 2 d ( . , n d^\ 

^Teae [ sme w ^9—39 sm d m 



1 



sin 2 (9 



+ (p 2 - 2)£M 2 $ = 5$. (46) 



Introducing the "tortoise" coordinate r* , which is de- 
fined by the relation dr* = -jjpj, together with a new 

wave function $ = *, one can reduce the wave equa- 
tion to the Schrodinger wavelike form with an effective 
potential V(r, S), 



where 



d 2 <$> 

dr* 2 



V(r,S) 



+ (uj 2 -V(r,S))<Z> = 0, 



2M 



2M 



(47) 



(48) 



The latter effective potential is identical to that of the 
Schwarzschild, up to the different angular eigenvalues S, 
which will be calculated in the next section. 
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FIG. 10: Effective potentials for a) I = 1, m = 0, S = 0, 1/10, 2/10 (from the top to bottom), b) £ = 1, m = 1, £ = 0, 1/10, 2/10 
(from bottom to top). The potential barrier is raised by the tidal force for non-zero m = ±1 and is lowered at the lowest modes 
(m = 0). 



VI. ANGULAR WAVE EQUATION 



The eigenvalues of the angular equation (|46|) can be 
calculated via two methods: the convergent and accurate 
Forbenius method and the method of expansion into the 
Associated Legendre Polynomials. Here we shall imply 
that the maximum of the effective potential is situated 
at its Schwarzschild distance r max given by (|42j) for each 
value of the multipole number £, i.e. we shall neglect 
tiny displacement of the peak induced by the tidal force 
£ . This approximation is fully justified as the inclusion of 
{£ , B) corrections to the maximum's location contributes 
only at the subdominant orders, as discussed in Sec. VII. 



Frobenius Method 



Equation (|46|) is reduced to the Heun's equation by 
setting x = sin 2 9 and thus giving 



transform equation (I49|) into the Heun's equation 



dx 2 



7 



x x — 1 x — an J dx 



d^ 



a/3x — q 



with 



x(x — 1) (x — an) 



$ = 0, 



(50) 



7 = m + 1, 5 = -, k = 1, 
m 2 +m + S + 2m 2 a + 2ma 



2 ' 



= 



4(7 

m + 3 



Off = 



and 



S'a; - 



w = 7 + £ — 1 = a + /3 — k = m + — . 



d 2 <I> 2 - 3x - 4a;cr + 5:r 2 cr $ n. 

2(x — jWra — 1) ~dx~ 4(x — l)x 2 (ax — 1)^ = ^ e solution to the Heun's equation can be expanded as 



where 



(49) 



follows 



+oc 



a = (p 2 - 2)£M 2 



ip(x) = ^2 c n y n (x), 



The above expression for a will be considered as a con- 
stant in the further calculations, which is justified by the 
approximation (|42p . In this way, we are only able to find 
even eigenfunctions for 8 — > —9. The odd eigenfunctions 
must be found separately. Setting <& = x^~ip( x ): one can 



if c n satisfies the recurrence relations 

F n c n -\ + G n c n + H n c n+ i = 0, 

where 



(n + a- l)(n + j3 - l)(n + 7 - l)(n + ui - 1) 
(2n + w-2)(2n + w-l) ' 



«n(n + w)(7 — 5) + (n(n + w) + aj3)(2n(n + lj) + j(oj - 1)) . . 
G„ = — — : , w „ — : — h n.(n + w )a ff + 5, 



#n = 



(2n + w + l)(2n + cj - 1) 
(n + l)(n + w — a + 

(2n + w + 2)(2n + cj + 1) : 
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with n = ^j 11 . 

Since the coefficients of c n -i, c„ and c n +i vanish indepen- 
dently, the eigenvalue can be obtained by taking G„ = 0, 

£(1 + I) (3 - U - M 2 + 2£a + 2£ 2 a - 6m 2 a) 

(-l + 2*)(3 + 2*) ■ 

(51) 

The above formula has been derived here for the " even" 
modes, i.e. modes with even values of I — \m\, though, 
using the expansion in terms of the Associated Legendre 
Polynomial, we show that this formula is valid also for 
odd eigenfunctions. 



force in the potential maximum. This means that instead 
of Eq. ([42)1 for p, one could use a general solution for the 
maximum of the effective potential (|48|l : 

p = -3 + 35 + v/(9 + 5(14 + 95))/25 (55) 

Comparison of the QNMs computed with the 
Schwarzschild values of p and with (|55|) for lower 
values of I, reveal a very small difference of about 0.01 
percent, which cannot be distinguished in figs. (fTTjl and 
(fT2")l . In the next section, the more accurate value given 
in (|55|) shall be used. 



B. Expansion in terms of the Associated Legendre 
Polynomial 

Introducing a new variable i in a different way, x — 
cos 9, equation ([4*o]) can be written as 



d 2 3> 
dx 2 



2x (2x 2 a - 2a + l) 



(x - l)(x + 1) (x 2 o - a + 1) dx 



-m 2 x 2 a + Sx 2 



(x - l) 2 (x + l) 2 (x 2 a - a + 1) 



-$ = 0. 



(52) 



Since the Preston-Poisson solution is accurate through 
the order of £ , one can expand the angular eigenfunction 
function and the corresponding eigenvalues in terms of 
a = (p 2 - 2)EM 2 as 



\?>\m\ 

5 = £{£+!) + crS', 



(53) 



(54) 



where £ = 0,1, 2, and m = -£, -I + 1, , 1-1,1. 

Substituting equations (|53|) and ([54)) into the wave 
equation (|52j) . one can obtain the 5' values, which are 
summarized for lower I and m in Table IPV1 



< = 





1 


2 


m — 





-10 1 


-2-10 12 


S' = 





4 8 4 

5 5 5 


24 12 24 12 24 
7 7 7 7 7 



TABLE IV: The eigenvalues of I = 0,1 and 2. 

Table HVl shows that the formula (|51l) is indeed valid for 
the low laying even and odd modes. With this data for 
the angular eigenvalues at hand, we can start numerical 
calculations of the quasinormal modes. 

In principle, one can have a slightly better approxima- 
tion of the estimation of the angular eigenvalue 5, if the 
tiny shifts of the value of the radial coordinate r in the 
maximum of the effective potential due to nonzero values 
of £ are taken into account, when calculating the tidal 



VII. ESTIMATIONS OF THE QUASINORMAL 
MODES 



An essential moment in the determination of the quasi- 
normal modes is the agreed notion of asymptotical re- 
gion, " infinity" , which is simply a distance much larger 
than the black hole radius. Once the near horizon and far 
regions are defined, we can formulate the boundary con- 
ditions for quasinormal modes. Quasinormal modes are 
solutions to the wavelike equation (|4T[) which are purely 
ingoing waves at the event horizon and purely outgoing 
waves at the far asymptotic region: 

$ 



far asymptotic region (56) 



Thereby, no incoming waves are allowed either from the 
horizon or from the far region. This means that quasi- 
normal modes correspond to the proper oscillations of the 
black hole's response to the external perturbation at late 
time, that is, when the source of the initial perturbation 
does not act anymore. In other words, the perturbation 
is considered as a " momentary" . 

In order to calculate quasinormal modes defined above, 
one can use the WKB formula, which takes into consider- 
ation the behavior of the perturbation near the top of the 
effective potential. The asymptotic region is located at a 
distance, which is at least one order larger than the black 
hole radius, such that one could qualitatively model the 
processes far from the black hole. At the same time, such 
a defined "infinity" must be located at a distance which 
is quite a few times smaller than the radius of the giant 
torus. The latter is necessary to avoid influence of the 
boundary effects near the torus of matter. Thus, there 
should be a "valley" far from the black hole and far from 
the edge of the torus, r^ -c r <C a which could model an 
asymptotic region. 

The QN spectrum is one of the essential characteristics 
of a black hole because it depends only on the black hole 
parameters and not on the way by which modes were 
excited. Supposing "J ~ e~ luJt , quasinormal modes can 
be written in the form 

U! = LOR e — UJj m i, 
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where positive wj m is proportional to the decay rate of a 
damped QN mode. The low laying QN frequencies have 
the smallest decay rates in the spectrum and thus dom- 
inate in a signal at a sufficiently late time. They can be 
calculated by the WKB approach [Hj], [ll]- Introducing 



Q 



V, the 6-th order WKB formula reads 



*<9c 



i=6 



i=2 



71 = 0,1,2. 



(57) 



where the correction terms were obtained in [24( , [25J ■ 
Here Q l means the i-th derivative of Q at its maximum 
with respect to the tortoise coordinate r* , and n labels 
the overtones. The WKB formula (|57|) was effectively 
used in a lot of papers (see [26| and references therein). 

As an additional check of the WKB data, we used the 
time domain integration scheme, which was related, for 
instance, in [27| and used in a number of other works (see 
[28| and references therein). The obtained time domain 
data shows about a 0.1 percent difference with the WKB 
values in the worst cases when £ = 1 and £ is moderate. 
The usage of the WKB formula is quite efficient here, 
because it has a very good accuracy for £ > n, while 
the potentially "unsafe" WKB mode £ = n = is not 
effected by the tidal force in the used approximation and 
coincides with the Schwarzschild one. 

The low laying quasinormal modes for various values 
of I and m are shown in figs. (fTTj) and (TT2"j) . The modes 
with I = to are noticeably different from those with £ > 
to: while the first are monotonically increasing when £ 



increases, the second one decreases. Such behavior of 
the spectrum can be explained as follows: the height of 
the potential barrier is lowered by the tidal force for £ > 
to, which means that it becomes easier penetrable for 
waves, leading to smaller real oscillation frequencies and 
longer living modes. At the same time £ = to modes 
are governed by the potentials raised by the tidal force, 
so that both real and imaginary parts of u grows, as £ 
increases. The lowered potential for all modes, except the 
one with the highest azimuthal number m = £, can be 
explained by a smaller resultant gravitational attraction 
near the top of the potential due to the tidal gravity, 
which is opposite to the gravitational attraction of the 
black hole. The non-monotonic " splash" of the imaginary 
part of u> in fig. (| 1 1 1) could probably be explained by 
approaching the limit of validity £ <C 1. 

In the limit of large multipole numbers f » m » 1, 
which is the eikonal regime, the peak of the effective po- 
tential is situated at 

r - = 3M - 3(2-7A^^ + ° (r4) - (58) 
Then, the first order WKB formula, 



(59) 



expanded in powers of l/£, gives: 



1 - l£M 2 
3V3M 



- ) —i \ n + - 
2 \ 2 



•O(l/0- (60) 
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In the limit £ = 0, the above expression reduces to the 
well-known form for the Schwarzschild black hole [29j . 
The quasinormal modes do not depend on m in the first 
two dominant orders: the dependence on m is in terms 
of order ~ 1/1. The imaginary part of ui does not depend 
on £ in the eikonal regime, which is in concordance with 
the very slight dependence of the imaginary part on £ 
in figures (2) and (3). It is well known that the eikonal 
formulas similar to (|60|) work very well already for mod- 
erate values of I and, in the eikonal regime, a universal 
behavior for all bosonic fields takes place. Therefore, ex- 
pression (p0|) is likely to be valid not only for scalar but 
also for other boson fields in the eikonal regime. 

An interesting question is the gravitational stability of 
the considered system of a torus and a black hole. The 
full assessment of the stability issue should be performed 
by considering the gravitational perturbations of the sys- 
tem, yet, as we can see from the above data, the scalar 
field shows no unstable modes. A tendency to the in- 
stability shows I = m modes through the decreasing of 
the damping rate, which however cannot be extended to 
sufficiently large values of £ as our effective potential was 
obtained for the regime of relatively small tidal force. 



VIII. CONCLUSIONS 

The black hole deformed by the surrounding matter 
and the magnetic field is a background with a relatively 
low symmetry, so that one could not expect the sepa- 
rability of variables in either Hamilton- Jacobi or Klein- 
Gordon equations. Using a few justified approximations, 
we have decoupled variables in both equations and made 
the problems of particle motion and quasinormal modes 
solvable with relatively simple analytical and numerical 
tools. 

Initially we have studied motion of charged, massive 
particles in the equatorial plane. Two approximations 
were used here: 1) we considered only a limited spacial 
region which starts at the event horizon and finishes at 
some large distance from the black hole, being still quite 
far from the place where the tidal force becomes signif- 
icant; 2) we were limited by relatively small values of 
the magnetic field and tidal force. These approximations 
were stipulated by applicability of the perturbatively ob- 
tained Preston-Poisson metric Q and was, therefore, un- 
avoidable. Consideration of motion only in the equato- 
rial plane was due to the inseparability of the Hamilton- 
Jacobi equation in the general case. This limitation is 
also justified, because, as it was shown in [3] for the Ernst 
solution, motion in the equatorial plane would be stable 
against small perturbations perpendicular to the plane, 
so that the particle which started to move in the equa- 
torial plane, would continue moving in it. However, an 
analysis of the stability against non-equatorial perturba- 
tions, though similar to [3j, is lacking. 

Motion of particles is qualitatively different for left and 
right handed rotations due to the opposite direction of 



the Lorentz force: when the signs of the particle's charge 
and angular momenta coincide (Larmor motion) we have 
a kind of cyclotron rotation in the magnetic field "per- 
turbed" by a black hole, while in opposite rotation (anti- 
Larmor motion), the presence of the black hole is essen- 
tial. We have found the energy and angular momentum, 
the binding energy and the region of radial stability on 
circular orbits for the Preston-Poisson space-time. From 
this we conclude that the tidal force, as well as the mag- 
netic field, considerably enhances the release of the bind- 
ing energy and makes the region of stability of circular 
orbits closer to the black hole. 

The second part is devoted to perturbations of the 
massless scalar field in the Preston-Poisson background, 
namely, to the estimations of the characteristic quasinor- 
mal modes. As the decoupling of variables is impossible 
for the Klein-Gordon equation in the general Preston- 
Poisson space-time, we used an additional approximation 
based on the fact that the low-laying quasinormal modes 
are "localized" near the peak of the effective potential. 
Therefore, we performed a kind of averaging of the tidal 
force by its value at the peak of potential. We found a 
significant decrease in the real oscillation frequency with 
the tidal force. A simple analytical expression for the 
frequency has been obtained in the eikonal regime. Tak- 
ing into consideration possible similarities between the 
quasinormal spectrum of a scalar field with the spectra of 
gravitational and other long-range neutral boson fields, 
we neglected the influence of the magnetic field when 
studying the quasinormal modes. Although we managed 
to estimate the modes with £ > 1, the monopole mode is, 
unfortunately, undistinguishable from its Schwarzschild 
value in the approximation. All the limitations of our 
analysis are a price for the simplicity afforded by the 
approximations employed. A more accurate numerical 
investigation would require very time consuming compu- 
tations and mathematical modeling of the above system. 
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